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1 , Problems 

No problems are .impetling progress, 

2 . Approach 

Two now approaches were developed during the preceding 
quarter (1 Jun,-3l March H'S!) and desesribed in Progress Report 
(to which referenice should be made in evaluating this report) . 
We now describe another now approach that has been devol-' 
oped during the current report period. The theoretical and 
practical importance of the "unsigned magnetic flux" or pole- 
strength of the earth, defined at radius r and time t by 

211 It 

P(r,t) J |* 1 (r, 0 ,t) I r^sinO dOdtli (1) 

*0 ‘0 

stems from two considerations. Firstly, this somewhat bi/sarre 
magnetic quantity is unique in that it is rigorously invariant 
in time at the top of earth's liquid core provided the core 
moves like a perfect conductor, which it must do on suffici- 
ently short time scale, since ohmic diffusion is not 

instantaneous. Secondly, the constancy of P at r ~ r « 

o 

3485 km provides one constraint for the magnetic forecast 
problem, since a mainfield model for some future epoch should 
produce the same value of P at r - r^. 

Because of these faces we believe it important to under- 
stand P both physically and mathematically. We have therefore 
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developed a method for evaluating P approximately at the surface 
of the earth r - a = G371.2 km (where P is also of interest 
since even though it is not constant, its time dependence is 
weak and smooth — see Figure 1 in section 3 below) . 

Expressing B as the curl of a vector potential A and 
then applying stokes' tiheorem shows that an alternative expres- 
sion for P is 

P(r,t) « 2 A(r ,0,(j),t) • d£ , (2) 

C 

where the contour integral is around each curve C on the 
spherical surface of radius r where the vortical magnetic 
field, B., vanishes. At the surface of the earth there is 
only one such curve, the magnetic equator, and it is every- 
where no more than about 13° displaced from the geographic 
equator. Consequently, small angle approximations to trigo- 
nometirc functions can be introduced and these lead to a very 
great simplification of the required integration. The details 
would be inappropriabe to give hero, but the first order re- 
sult is simply that 

P(a,t) « -4'ira^(gj^^ - 

vdiere g^^^, axial dipole and axial octupole Gauss 

coefficients. The above approximation is found to be accurate 
to about 0.6% by comparing its prediction with an accurate 
numerical calculation. A better approximation is being 
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developed, but it involves a fairly heavy amount of tedius 
algebra . 

A different but simple result of beauty is obtained 
for the conceptually important case where the magnetic field 
is purely a centered (but tilted) dipole. Then in tilted 
dipole coordinates the value for (~ dipolar contribution 
to P) can be worked out with ease, the result being 

3 

= 4n ^ _ ,3, 

Stated in words: the unsigned dipole magnetic flux is an 

inverse function of radius proportional to the magnetic 
dipole moment. At earth's surface, r = a, the unsigned magne^ 
tic dipole flux per unit area is the magnetic dipole moment 
per unit volume. Equation (3) reveals that P^, and therefore 
P, depends non-linear ly on each Gauss coefficient in general. 
This conflicts with some recent work of Hide and Malin (1981, 
Proc. Roy. Soc. A374 , 15-33) and may in part explain some 
mysterious time dependence in their calculations. 

3. Accomplishments 

We focus on two particular accomplishments, both involving 
the pole-strength as discussed above. 

a ) Magnetic determination of the depth to the core - 
mantle boundary using MAGSAT data . 

Since P generally changes in time except at the top of 

earth's core, r = r_, a calculation of P as a function of 
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depth £or two different times, t^, t^, provides a magnetic 
determination of the radius of the core-mantle boundary. 

This idea was proposed by Hide and applied by Hide and Malin, 
but they used a mainfield model together with a s(?.eular varia- 
tion model and obtained results which are suggestive but not 
definitive . 


We have produced what we believe to be a definitive 
test of the concept by reducing the dependence on secular 
variation (which is highly uncertain when evaluated at depths 
within the oarth) to nec,ligib.lo proportions. For t^ we 
choose 1965 and adopt the so-called definitive mainfield model 
of Barrae lough, Harwood, Leaton, and Malin. It is an intor- 
polative rather than extrapolative model built on a i^assive 
data set centered in time near 1965. Most data used wore 
within 5 years of this epoch so minimum reliance on a SV 
model was required in bringing the data to conm\on epoch. 

For t 2 we choose 1900 and adopt the mainfiold model from 
MAGSAT of Langol ot al. MGST 6/80. This used data only from 
the two quiet days 5,6 Nov. 1980, so it too is virtually free 
of secular variation. This model was a fit to the data at 
N == 13, but we have truncated back to N » 8 since that was 
the value used for the 1965 model (we plan to re-do the 
calculation using an N ~ 8 fit to the MAGSAT data, but expect 
the results to change negligibly) . 

The primary result is plotted as Figure A. It shows 
how P (r ,t 2 ) “P (r ,t^) varies with depth in the mantle (assumed 
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wi t h 1 n 1 . n 3 V. o 1 1 1 s accurat e 
is hiqhly atqiuflcant in tnir viinv. 


b ) H i a t iir um 1 4 y\f t jie uns ii^ju?d 

nui^nel ic I'Uik^ Ih 

As planned in thv' K\st proqri'ss refurirt , we Inive now 
completed numerical v'valiu\l‘ ion oi P as t\ 1 unction of time f(vr 
r ' a t37l,3 km. At earth’s surface, the hiqher liarmonics 
ol the aeomaqnel ie field do nob eontribute stronqly to the 
valiu.^ ol P (as they do at earth's eoi*e) . Therefore, it is 
meaninqful to use hisborical models of the aeomagnetic field 
from times back a century or more ago, for which only the 
fairly lew order harmonics a.re available. 

Figure 1 shows the time dependence of P(a,t) using the 
N - 4 models constructed by Barraelough (1974) for 50 year 
inter\Mls back t.o ItOO, One sbandard deviation error bars 
are also included. Kvidently for about the last 200 years, 
P(a,t) has been declining steadily, smoothly and novvrly 
linearly. It would make sense to predieb that at 1985, 

P(a,b) will be about 1.5 8 GWb (for N 4). (The dot on this 
fiaure at M 8 shows the value calculated from MGST 6/80, 



FIEUREi 1: 
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thoroby indicating the relatively small contribution of the 


higher order harmonics) . 


• Signifi ca nt Resul ts 

MAGSAT data has been used, in conjunction with a high 
quality mainfield model for epoch 1965, to confirm that the 
radius of earth's corc-mantle boundary can indeed the found 
magnetically with impressive accuracy. The unsigned flux 
linking earth's core and mantle should therefore be a legiti- 
mate invariant for a span of time. The value from MAGSAT of 
this constant is 16.056 GWb (gigiwebers) . 

• i'ub i Ic a ■ 

(Related to this project, but not directly supported 
by it.) 

Mdward R. Renton, "A simple model for determining the vertical 
growth rate of vertical motion at the top of earth's 
outer core," Physics of the Earth and Planetary Interiors 
24 (1981) , 242-244. 

Lor ant A. Muth and Edward R. Benton, "On the frozen flux 

velocity field at the surface of earth's core necessary 
to account for the poloidal main magnetic field and its 
secular variation," Physics of the Earth and Planetary 
Interiors ^ (1981), 245-252. 
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6 . Recommendations 

None . 

7f Funds Expended Through May 31^ 1981 ; 

8. Data Utility 

It should be apparent that the MAGSAT data as acquired 
and modelled is of very great utility to this project. 


